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We study theoretically the entanglement of two-photon states in the ground state of the inter- 
subband cavity system, the so-called polariton vacuum. The system consists of a sequence of doped 
quantum wells located inside a microcavity and the photons can interact with intersubband exci- 
tations inside the quantum wells. Using an explicit solution for the ground state of the system, 
operated in the ultrastrong coupling regime, a post-selection is introduced, where only certain two- 
photon states are considered and analyzed for mode entanglement. We find that a fast quench of the 
coupling creates entangled photons and that the degree of entanglement depends on the absolute 
values of the in-plane wave vectors of the photons. Maximally entangled states can be generated by 
choosing the appropriate modes in the post-selection. 



I. INTRODUCTION 

With the advent of quantum information theorjP, the 
phenomenon of entanglement not onl y re mained a mys- 
terious feature of quantum mechanics^!, but became a 
resource to perform tasks that are not feasible with clas- 
sical resources. Examples are quantum communication 
protocols, which make use of entangled states like quan- 
tum key distribution^, quantum teleportatiorP or super- 
dense codinjP, or the realization of a quantum repeated. 

Entangled photon states are often used to implement 
the protocols mentioned above. Today, there exist several 
different proposals for the production of bipartite entan- 
gled photon states, most prominently type-11 paramet- 
ric down-conversiorP and biexciton decay in a quantum 
doP. 

The fundamental requirements for such a photon-pair 
source to be used in quantum information processing are 
that the states have to possess a sufficient amount of 
entanglement and that the production of the two photons 
has to be deterministic and efficient. Determinism means 
that the release of the photons can be triggered by some 
external control parameter. Efficiency means that the 
probability for this event is near unity. 

Here, we study the intersubband cavity system, for 
which the emission of correlated photon pairs was pre- 
dicted theoreticallylSl and can be triggered by modu- 
lating the light-matter interaction between microcavity 
photons and electronic excitations in the quantum wells. 
Those intersubband transitions are mainly used in quan- 
tum well infrared photodetector 3^ and quantum cascade 
lasers^tt^. Embedded in a microcavity, it is possible to 
reach a regime of ultrastrong light-matter coupling^HIOD 
in which the vacuum-field Rabi frequency can be of the 
order of the intersubband transition frequency and the 
ground state of the system, a squeezed vacuum, con- 
tains already a non-zero number of photons. Another 
type of system, which can reach the ultrastr ong co upling 
regime as well, are superconducting circuits22H22l ) where 
the emission of quantum vacuum radiation was just re- 
cently demonstratecpH 

In this paper, we analyze the ground state of the in- 
tersubband cavity system, the so-called polariton vac- 
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FIG. 1. (a) The system under consideration: A sequence of 
doped QWs (gray) inside a microcavity of length L c (cav- 
ity mirrors in blue) form the intersubband cavity system. 
The light-matter interaction in such a system depends on the 
propagation angle 6 of the cavity photons. The transverse- 
magnetic polarization is indicated by the electric and mag- 
netic field vectors. The electric field E lies in the plane of 
incidence and the magnetic field B is perpendicular to it, 
and both are perpendicular to the photon wave vector fc to t- 
(b) The plane of incidence to better demonstrate the TM- 
polarization of the photons. 



uum, related to two-photon entanglement. We use an 
explicit expression for the polariton vacuum and, after 
post-selecting certain photonic states, quantify the mode 
entanglement between the photon pairs via the concur- 
rence. 



II. THE INTERSUBBAND CAVITY SYSTEM 

The intersubband cavity system was studied theoret- 
ically by Ciuti et aZP^. It consists of tiqw identical 
quantum wells embedded inside a semiconductor opti- 
cal microcavity. The quantum wells are assumed to be 
negatively charged with a two-dimensional electron gas 
(2DEG) with density A^deg that populates the first sub- 
band. We consider the interaction of intersubband exci- 
tations between the two lowest subbands and photons of 
the fundamental cavity mode. 

The ultrastrong coupling regime, in which the vacuum- 
field Rabi splitting is of the order of the intersubband 
transition energy, can be reached due to the large dipole 
moment^ of intersubband transitions and the collective 
coupling to all electrons of the 2DEGs. In this regime, the 
rotating wave approximation^ is not valid anymore and 
the full light-matter interaction Hamiltonian including 
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FIG. 2. Subband energy structure of a quantum well (QW) 
formed in a semiconductor heterostructure. (a) In real space, 
along the growth direction z of the structure, the semicon- 
ductors forming the QW are denoted as A and B. The QW 
contains a two-dimensional electron gas (2DEG). Here, we 
study intersubband transitions between the first two subbands 
n = 1 and n = 2 with transition energy huj^. (b) The same 
situation in k space, where the 2DEG populates states up to 
the Fermi energy Ef, with a corresponding wave vector &f- 



the antiresonant terms has the forrrP^ 



H = H + H TC 



which consists of the three terms 



(1) 



H Q = ^c(fc) [a{a k + 2 J + E ^™ h lKi ( 2 ) 
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+ hD ( k ) ( a l a k + a k a l j 
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(3) 



(4) 



where Hq describes the uncoupled photonic and elec- 
tronic systems, H rcs is the resonant part of the light- 
matter interaction and the antiresonant terms, usually 
neglected under the rotating wave approximation, are 
given by £f a nti- 

The operator a k ^ annihiltes (creates) a cavity photon 
with in-plane wave vector k = (k x ,k y ) and transverse- 
magnetic (TM) polarization. The reason why only this 
photon polarization is considered, is the selection rule for 
intersubband transitions 1 -^. The dipole moment points 
along the growth direction, so the exciting radiation must 
have a finite electric field component in z-direction. As 
can be seen from Fig. [T] TM-polarized light where the 
magnetic field component is perpendicular to the plane of 
incidence, has an electric field component in z-direction 
when it propagates under a finite angle 9. However, if 
the polarization is not purely transverse-magnetic, it is 
still sufficient to consider only one polarization direction, 



since only the TM-polarized part of the radiation is ab- 
sorbed. 

The only bright collective intersubband excitation with 
in-plane wave vector k is described by the operators b k 
that fulfill Bose commutation relations — Sk.k' 

in the weak excitation regime, i.e. when the number of 
intersubband excitations is much less than the number of 
electrons forming the 2DEGP^. 

The dispersion of the fundamental cavity mode ui c (k) 
is given by 



(5) 



where c is the speed of light, e is the dielectric constant of 
the material used as cavity spacer and the quantization 
of k z can depend in a complicated way on the boundary 
conditions. In the following, k = \k\ is the length of the 
in-plane wave vector. 

The vacuum Rabi frequenc y £l R (k) for the intersub- 
band cavity system is given b} * 10 * 26 ! 

i 

e 2 iV 2D EG^Qw/l2 a; 12 



Q R (k) 



2e em Lfuj c (k) 



sin 2 6(k) 



(6) 



Here e is the elementary charge, So is the vacuum per- 
mittivity and mo is the free electron mass. L° denotes 
an effective cavity thickness that depends on the type of 
cavity mirrors and Wq^, is an effective number of em- 
bedded quantum wells since not all quantum wells are 
equally coupled to the photon field. f 12 is the oscil- 
lator strength of the transition between the subbands. 
For a deep rectangular well, one can use the approxi- 
mation /12 ~ mo/m* (m*\ effective electron massjM^M 
Finally, 0(k) is the propagation angle of a cavity photon: 
sin#(fc) = kj ' \Jk 2 + fc 2 . The dispersive coupling param- 
eter D(k) can be approximated by D(k) ~ f2 2 l (fc)/w 12 , 
which is valid for deep rectangular wells and exact for 
parabolic well potentials^. 

The Hamiltonian ([I]) can be diagonalized with an ex- 
tended Bogoliubov transformation^, also known as Hop- 
field transformation^, where new bosonic operators 

P 3 ,k = Wj (k)a k + Xj (k)b k + y {k)a}_ k + z {k)W_ k (7) 

are introduced that describe a quasiparticle called inter- 
subband cavity polaritoii 2 ^! and j indicates wether it be- 
longs to the lower (j = LP) or upper (j — UP) polari- 
ton branch. The wave vectors k are still meant to be 
in-plane. By an appropriate choice of the Hopficld coef- 
ficients Wj(k), Xj(k), Uj(k) and Zj(k), which are already 
taken to depend only on fc, the Hamiltonian becomes di- 
agonal 



H = Er 



E E 

je{LP,UP} k 



(8) 



Here, Eq denotes the ground state energy. The result- 
ing lower and and upper polariton dispersion are shown 
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FIG. 3. The polariton energy dispersions hLJLp(k) and 
hu>\jp(k) as a function of the absolute value of the in- plane 
wave vactor k. Here, ^qw = 

50 GaAs/AlGaAs (e = 10.0, 
/12 = 14.9) QWs are assumed to be located inside a cavity of 
length = 2 fjxa and doped with a two-dimensional electron 
gas density of A^deg = 10 12 cm -2 . The lower plot shows the 
anticrossing region where the dotted lines are the bare cav- 
ity photon dispersion hu) c (k) and the constant intersubband 
transition energy huii2 = 150 meV. The vacuum-field Rabi 
splitting is 2h£lB_(k Ica ) . 



in Fig. [3] Here, as well as for all further quantita- 
tive resultes, we assume GaAs/AlGaAs quantu m wells, 
which have been commonly used experimentall y 16 * 29 " 31 !. 
Hence, the material-dependent parameters are /12 = 14.9 
( TO GaAs = 0-067) and e = 10. Furthermore, the number 
of embedded quantum wells tiqw, the length of the mi- 
crocavity L c and the intersubband transition frequency 
W12 , which is determined by the quantum well depth and 
thickness, can be adjusted during the manufacturing pro- 
cess. The density of the two-dimensional electron gas 
can be varied experimentally. To obtain the results of 
Fig.|Z| we chose n^ w = 50, L° ff = 2 fim, hui 12 = 150 me V 
and iV 2 DEG = 10 12 cm~ 2 as one particular set of experi- 
mentally reasonable values of the parameters mentioned 
above. 

The ground state of the intersubband cavity system, 
operating in the ultrastrong coupling regime, is not the 
ordinary vacuum |0) , in which there are no cavity photons 
and no intersubband excitations present 



o fc |0) = h fc |0) = 0, 



(9) 



but a state |G) that exhibits no intersubband cavity po- 
laritons 



p j!k \G)=0, je{LP,UP}. 



(10) 



Without knowing the explicit form of | G) , one can show 
that the ground state has some peculiar properties in 
the ultrastrong coupling regime that were worked out in 



Ref. [TUl whereof the essential ones are that it contains a 
finite number of photons: 

<G|a t fe a fc |G) = |y L p(fc)| 2 + |j/up(fc)| 2 , (11) 

and photons with opposite in-plane wave vectors k and 
— fc are correlated: 

(G\a k a^ k \G) = -w£p(fc)y LP (fc) - w{jp(fc)j/ UP (fc). (12) 

One can see that only if the light-matter interaction is so 
strong that the Hopfield coefficients yi,p(k) and y\jp(k) 
are reasonably large, \yj(k)\ 2 ~ 0.1 (i.e. when the antires- 
onant terms Q cannot be neglected and therefore the ex- 
tended Bogoliubov transformation (7| is necessary), the 
ground state \G) differs significantly from the vacuum 
state |0). 



The idea is now that the correlations ( 12 ) can lead to 
entanglement of two photons propagating in opposite di- 
rections. These phot ons are, however, virtual excitations, 
but it is conjecture d 10 * 32 * 3 - 3 -* that they can be released 
by a non-adiabatic switch-off (quench) of the vacuum 
Rabi frequency f2;a(fc). An experimental approach to this 
scenario is an ultrafast cha nge of th e density A^deg of 
two-dimensional electron ga d 16 * 30 * 3 -U. One mechanism to 
achieve a modulation of the parameter A^deg is a gate 
voltage, which can lead to the depletion of the QWiP^. 
The rapidity is restricted by capacitance of the gates, 
however. Another implementation uses two asymmetri- 
cally coupled QWs, in which one QW can be charged 
by electron tunneling and this process can happen on 
the picosecond time scale or faster^. A promising idea 
to achieve an ultrafast coupling modulation is an all- 
optical control scheme^, in which electrons from the va- 
lence band are resonantly excited to the first subband 
by a femtosecond laser pulse. In this manner, it could 
be demonstrated experimentally that the coupling be- 
tween the cavity photon field and the intersubband tran- 
sitions in the quantum wells can be switched on in a 
time shorter than a cycle of light in the microcavity. In 
Ref. |3H and [33J the spectrum of the radiation exiting the 
cavity was derived in more detailed calculations, when a 
time-dependent coupling i7p, (fc,i) is predominant in the 
system and it is predicted that the vacuum radiation rises 
above the black body radiation. 



III. EXACT GROUND STATE 

A pioneering calculation of the polariton ground state 
of a bulk dielectric was given by Quattropani et aZPSl. 
The solution is given by independent photon and polar- 
ization states. Since the Hamiltonian of the intersubband 
cavity system is similar to the one in Ref. I28| we use their 
treatment to determine the explicit form of the polariton 
vacuum |G) being the ground state of the Hamiltonian 
|TJ of the intersubband cavity system. The difference is 
just that the sums in ([I]) cover all in-plane wave vectors 
and the intersubband transition frequency w 12 is taken 
to be dispersionless. 
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The ansatz for the polariton vacuum \G) is: 

1 hJ2[G(k)(aial k +blbl k )+F(k)(albl k +b{al k )} 
\G) = ^ e k |0)- 

(13) 

N is a normalization constant and the expansion coeffi- 
cients G(k) and F(k) have to be determined in order to 
satisfy the definition of the polariton vacuum (10): 



Pj,k\G) = 0, je{LP,UP}. 

We anticipate that the functions G(k) and F(k) will only 
depend on the absolute value of the in-plane wave vec- 
tor k. After some algebra using commutation relations, 
which is explicitly given in Ref. 1281 the action of a k and 
b k on \G) is, e.g.: 

a k \G) = (G{k)al k + F(k)bl h ) \G), (14) 



b k \G) = (G(k)bl k + F(k)c 



\G). 



(15) 



Inserting (14) and (15) into the definitions ^ and (10), 
one obtains a system of equations for the coefficients G( k) 
and F(k): 



Wj{k)G(k) + Xj(k)F{k) + yj{k) = 0, 
Wj(k)F(k) + Xj {k)G{k) + Zj {k) = 0, 

which has the solutions: 

Xj(k)zj(k) — Wj{k)yj{k) 



G{k) 



w*(k)-x?(k) 



[1 w](k)-x?(k) 



(16) 
(17) 



(18) 
(19) 



That this can be fulfilled simultaneously by the Hopfield 
coefficients of the lower an d uppe r polariton, can be seen 
from the following relation^ 27 * 28 ! 

Whp(k) = xup(fe), x LP (k) = w V p(k), 
VLp(k) = zvp(k), 2Lp(fc) = yvp{k). (20) 

By inserting the explicit expressions of the Hopfield coef- 
ficients, the expansion coefficients can be rewritten as^ 

G , k , = g!2 + ^c(fc) - ^Lp(fc) - ^Up(fc) , 21 , 

W12 - (k) - w LP (k) - aj V p (k) ' 



(22) 



Finally, the polariton vacuum \G) is calculated to be 

1 iI2G(k)(ata ] .+b ] .b t 1 .-2i r ," 1 , 2 . - aW .) 
|^ _ _|_ g 2Z-- v )\ fc -fc^ k -h s2 R (fe) fc -hj ^23^) 

because the last two terms in the exponential can be 
combined, and with the normalization TV given by 



N = 



n(Kp(fc)i 2 +iz LP (fc)i 2 ) 5 . (24) 
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FIG. 4. The absolute values of the expansion coefficients G(k) 
and |F(fc)| as a function of the length k of the in-plane wave 



vector. (Parameter values: e = 10, /12 = 14.9, uq W 
Lf = 2fim, huj 12 = 150 meV, iV 2 DEG = 10 12 cm" 



50, 



The dependance of G{k) and |.F(A;)| on k is plotted in 
Fig. g One can see that |F(fe)| is about one order of 
magnitude larger than G(k) and 



G(fc)«|F(£:)| 2 «l, 



(25) 



in which this was checked for a wide range of experimen- 
tally acceptable values of the parameters riqw, L c , cj\2 
and TV2DEG- Therefore, the polariton vacuum state \G) 
will be expanded in a Taylor series to the second order 
in these coefficients as an approximation. 

If the light- matter interaction is turned off (ilpi(fc) = 
0), then G(k) and F(k) are zero, so the ground state 
would be the ordinary vacuum |0), as expected. 



IV. PHOTON ENTANGLEMENT 

As seen in the previous section, the intersubband cav- 
ity system contains a finite number of photons if it is 
in the ultrastrong coupling regime. One possibility of a 
photon pair in the ground state \G) to be entangled will 
be studied below. After specifying the t ype o f entangle- 
ment, it is quantified by the concurrenc e 3 * 35 l 



A. Mode entanglement 

We first define the type of photon entanglement. Since 
the transverse-magnetic polarization of the interacting 
photons is fixed by the selection rule for intersubband 
transitions^, polarization entanglement as achieved with 
parametric down-conversion or the biexciton decay is out 
of the question. But there exist anomalous correlations 



between photons with opposite in-plane momentum ( 12 ) 



(G\a k a- k \G) = -w£ P (fc)y L p(fc) " Wup(k)yvp(k). 

Our idea is to test the photonic states in \G) for mode or 
frequency entanglement!^. 
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FIG. 5. Two photons with different frequencies (colors) leav- 
ing the cavity in opposite directions. The two subsystems left 
(L) and right (R) are defined via the sign of k x . Since there is 
a difference in frequency, the photons have different in-plane 
wave vectors. 



In the following, we limit the treatment to only one 
direction, since the correlations ( 12 ) occur only for pho- 
tons with exactly opposite in-plane wave vectors. This 
direction is chosen to be the x-direction. This is shown 



schematically in Fig. [5 Photons with a negative (posi- 
tive) ^-component of the wave vector belong to subsys- 
tem L(R) for left (right). 

\G) itself is a vector from a Hilbert space T~L 1 which has 
a tensor product structure % — T a <g> T h , where T a and 
F b denote the Fock spaces of the photons and intersub- 
band excitations, respectively. The situation depicted in 
Fig. [5] is described by vectors in a subspace 'Hlr C H, 
which is itself a tensor product of "Hi, and Hr, the Hilbert 
spaces of the subsystems L and R, Wlr = Hi, ® Hr- 
Since \G) contains states from outside "Hlr, we project 
onto "Hlr with an appropriate projection operator, which 
will be described in the following and could be realized 
experimentally by a post-selective measurement. 



B. Post-selection 

The post-selection needs to fulfill the following require- 
ments: first, we only allow for states in which two pho- 
tons with opposite in-plane wave vectors appear. In addi- 
tion, the post-selection is even more restrictive in terms of 
allowed modes. We only consider two different modes in 
each subsystem L and R, respectively, that have, however 
the same absolute value of the in-plane wave vector. I.e. 
we choose k and q, with k =/= q, and consider the modes 
k = (k, 0) and —k = (— fc,0) and accordingly q — (q,0) 
and q = (—q, 0), where all the wave vectors point along 
the x-direction. So the basis states of "Hi, are 

|fc>L = a f _ fc |0) a , (26) 
|g>L=« t _ q |0) Q , (27) 

where |0) a is the photon vacuum. The basis of Hr is 



|fe>R 



at|0)„, 



oL|0).. 



(28) 
(29) 



Hence, one possible product basis of "Hlr is 

|fc) L ®|fe) R , |fc) L <g) |q) R , 

|q) L ®|fe) R , |g) L ® |g>R. (30) 

For all further calculations, the polariton vacuum \G) 
is expanded to the second order in the small expansion 
coefficient G(k), 



\G) 



1 

N< 
1 

N 
|G ( 



|0) 



2 ^ 



G{k)% 




2 )\ 



|0) 



(31) 



where N is a new normalization constant to preserve 
(G^\G^) = 1 and the operator Tj is 



71 



°fc°-fc 



2i 



W12 



(32) 



^R(fc) 

Thus, the two-photon states with opposite in-plane wave 
vector, i.e. the states that fulfill the post-selection re- 
quirements, in linear order are 



4 a -fel 



= |fe>L|fc>R, 

= \<l)h\q)R, 



(33) 
(34) 



and in second order 



blbl k ,\0) b , 



blbl^h 



"W- q bl 



(35) 
(36) 
(37) 
(38) 



|fc} L |fc>R 

I<?}l|<Z>r< 

|g>L|fc>R®&l fc 6ilo>i 

\kh\q)R®bl q bl\0) b 

Here, the explicit expression of the tensor product |0) = 
|0) o ® \0)b of the individual vacuum states for photons 
and intersubband excitations is used and k! can be an 
arbitrary in-plane wave vector. 

We carry out the post-selection by projecting onto 
these states with a projection operator "Plr, 

1 



LR) 



LR 



-Vlr\G( 



2 )\ 



(39) 



with Ni,R being a necessary normalization constant, since 
the operation is a projection. As an intermediate re- 
sult, we obtain the pure state IV'lr) in which all the 
two-photon states fulfilling the conditions of the post- 
selection are extracted. We give an explicit expression 
for |?/>lr) in Appendix A. 

As we will see below, the reduced density matrix g a 
of the photonic system is needed for the calculation of 
the entanglement. We compute g a by tracing out the 
intersubband excitations 



0° = Tr 6 |V>LR}^>LR| 

Z{k) Y(k,q) 

X(k,q) 

X(k,q) 

Y{k,q) Z(q) 



1 

Mr 



(40) 



G 



Here, the matrix representation is in the basis (301 and 
the abbreviations 



X(k,q) = \F(k)\ 2 \F(q)\ 2 , 
Y(k,q)=G(k)G(q)x 
1 



l+^-S)-\F(k)\ 2 -\F( q )\ 2 



Z(k) = X(k,k) + Y(k,k), 



(41) 

, (42) 
(43) 



were introduced and S is the sum over all expansion co- 
efficients squared 



(44) 



The value of S (see Appendix B) depends on the sample 
area A. In the limit of A ^> 4? — , one can take the limit 
S — > oo and obtains 



a 5— foo 



1 



G 2 {k) G{k)G{q) 




G 2 (k) + G 2 {q) I 
G(k)G(q) G 2 (q) 



This corresponds to the pure photon state 

m 1 



^/G 2 (k)+G 2 (q) 



(45) 

(G(k)\k) L \k) R + G(q)\q) L \q) R ). 

(46) 



C. Measure of entanglement 

The state g a , which we derived from \G^), describes 
two photons that propagate with opposite in-plane wave 
vectors in the microcavity and that can potentially be 
released by an appropriate time-modulation (quench) of 
the Rabi frequency il^ (k). Since one chooses the modes k 
and — k and accordingly q and —q via the post-selection, 
the photons effectively form a two-qubit system. For such 



a system, the entanglement for mixed states can be cal- 
culated analytically without evaluating a convex roof ex- 
plicitly from the density matrix by way of the concur- 
rence ( '( '. With this function, the so-called entangle- 
ment of formation Ep(gp^ of two qubits can be easily 
calculated vigP^ 

m.)->( !±^g@ ), (47) 

with the binary entropy 

h(x) = -xlog 2 (x) - (1 - i)log a (l - x). (48) 

The concurrence itself is given bySS 

C{q) = max{0, A x - A 2 - A 3 - A 4 }. (49) 

Here, Ai to A4 are, in decreasing order, the square roots of 
the eigenvalues of the matrix g£and q is a transformation 
of the density matrix given by^ 

Q = (<r„ ® a y )g>*{o y ® a y ), (50) 
where a y is the Pauli y matrix and the * denotes complex 
conjugation. 

1. Analytical results 



For the photonic state g a (40), the parameters Ai to 



A4 are found to be 



Xi = j^-(y/Z(k)Z(q) + Y(k,q)), (51) 

LR, 

X 2 = ^-(yZ(k)Z(q)-Y(k,q)), (52) 

(53) 



LR 

^3,4 = -^-X(k,q) 



N 2 



Hence from (49), we obtain in connection with (41 )-(43 ) 



for the concurrence 



C(g a ) = C(k,q) = 



N 2 ■ 



G(k)G(q) 



1+ S )-\F(k)\ 2 -\F(q)f 



-\F(k)\ 2 \F(q)\> 



(54) 



The concurrence thus only depends on the absolute val- 
ues k and q via the expansion coefficients, which were 



given in (|21| and (|22j), 

W12 + w c (k) - ^Lp(fe) - w UP (fc) 



G(k) 
\F(k)\ 



W12 - oj c (k) - u>Lp(k) - wup(fc) ' 
W12 



n R (*) 



G(k). 



To show the dependance of the concurrence on k and q, 
we have to evaluate the sum S explicitly. Using a sample 



I 

area A = (200 /1m) 2 , we find S = 0.716. The result is pre- 
sented in Fig. [6] for the same parameters as used before, 
where C(fc, q) is shown in a density plot as a function of 
the modes k and q. Below it, we show cuts for differ- 
ent values of q to illustrate the dependancy of the con- 
currence better. One can observe two branches of high 
entanglement that appear for large values of k and/or q. 
Their appearence can be explained by the characteristics 
of the expansion coefficient G(k) for large k. G(k) tends 
to zero as 1/fc 2 , hence |F(fc)| 2 scales as 1/fc 3 . Hence for 
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FIG. 6. (a) The concurrence C(k,q) as a function of k and 
q for GaAs/AlGaAs quantum wells. (Parameter values: e = 
10, fa = 14.9, n c cf w = 50, Lf = 2/xm, hu} 12 = 150 meV, 
JV2DEG = 10 12 cm -2 .) In (b), we plot C(k,q) for fixed values 
of q (given in the respective legend) and for a large range of k 
from to 3 ■ 10 7 m _1 , which corresponds to a photon energy 
of 1.9 eV (450 THz, red) and in (c) magnification of the range 
up to k = 0.5 • 10 7 m* 1 (320 meV, 80 THz, mid infrared). 



the diagonal branch, i.e. k q, the \F\ terms in (54 1 and 
( A2 ) can be neglected and we obtain 



C(k,q) 



2G(k)G(q) kxsq 
G 2 {k)+G 2 {q) ~ 



(55) 



Accordingly, photons in the visible regime are almost 
maximally entangled if their wave numbers are of the 
same size. 

The other branch appears if G(k) m G(q) and the 
modes are far from each other. We give a more precise 



analysis of expression |55j) in the next subsection, where 
the limit of large sample areas is worked out. 



2. Large-cavity limit 

The case of a large cavity, i.e. a large sample area A, 
is described by the limit S — > 00, see Appendix B. The 
concurrence in this case is calculated to be 



C(k,q) = 



2G(k)G(q) 
G 2 (k)+G 2 {q)' 



(56) 



We show the result in Fig. [7J As one can see, the concur- 
rence always has two maxima if q is held constant. One 
maximum appears for G(k) — G(q) and k 7^ q. There, 
photons are maximally entangled since we have C = 1. 
However, this maximum is relatively sharpely peaked 
and if one realizes the post-selection experimentally by 
choosing a certain finite fc-range, the entanglement will 
be reduced. The other maximum appears when k = q, 
which seems to be an artefact of the calculation, since this 
case was excluded in the calculations above. The reason 
for the exclusion is that the two-photon states would be 
seperable and hence not entangled. In this case, how- 
ever, the maximum is quite broad. Therefore, for a given 
q there exists a wide range of corresponding modes k, for 
which the two photons are almost maximally entangled. 
In the second plot of Fig. [JJ we show a magnification for 
wave vectors up to 0.5 • 10^n _1 to show the dependance 
of C(k,q) on k clearer around the first maximum, which 
is not visible in the previous plot. In particular at the 
intersubband resonance, which is around 0.2 • 10 7 m _1 , 
the two maxima approach each other so that by selecting 
different modes around the resonance, the entanglement 
of the photons can be made almost maximal. The corre- 
sponding photon energies are in the mid infrared, about 
150 meV. 



V. CONCLUSION 

An efficient and deterministic source of entangled pho- 
tons is needed in quantum information processing. In this 
work, we examined a new scheme of photon production, 
based on the emission of quantum vacuum radiation from 
the intersubband cavity system. Because the triggered 
photon emission is based on a non-adiabatic modulation 
of the system's ground state, an exact expression for this 
state could be used. Since the ground state consists of 
an infinite number of photonic and electronic states, we 
propose a post-selective measurement to reduce the pho- 
tonic system to an effective two-qubit system, in which 
the qubit state was defined as two different in-plane wave 
vectors. The so-called mode entanglement of the photons 
is quantified by the concurrence. We found an analyti- 
cal expression for the concurrence, which depends on the 
absolute values of the chosen wave vectors. We find that 
the concurrence, and therefore the entanglement of the 




0.2 0.3 
k[10 7 m" 1 ] 



of large sample areas, there exists a continuous set of 
mode pairs for which the concurrence is 1, i.e. the pho- 
tons are maximally entangled. Also, in this case it turns 
out that for photon energies around the intersubband 
resonance, which is in the mid infrared regime of the 
electromagnetic spectrum, the photons are almost maxi- 
mally entangled, the concurrence being close to 1. This 
is fundamentally important for the possible use in quan- 
tum information processing. Furthermore, a high degree 
of entanglement can be achieved if the modes chosen in 
the post-selection are close to each other. Therefore, one 
could extract entangled photon pairs in technologically 
relevant frequency domains like a telecom wavelength. 
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FIG. 7. The concurrence C(k, q) in the case of large sample 
areas for GaAs/AlGaAs quantum wells. (Parameter values: 
£ = 10, /12 = 14.9, Uq W = 50, Lf = 2/im, huj 12 = 150 meV, 
AT 2DEG = 10 12 cm- 2 .) 



post-selected photons, is non-zero. In the limiting case 



Appendix A: Post-selected state 

The pure state |^lr)j which contains all two-photon 
states fulfilling the conditions of the post-selection, is 
given as 



N L r 



LR 



(G(*)|fc> L |fc> R + G(q)\q) h \q) R ) ® ( |0) 6 + ± £ G(k')b{,bl k , \0) b j + J F 2 (fc)|fc) L |fc) R ® bl k b{\0) t 



+ F 2 (q)\q) L \q)K <8 bl q b\\0) b + F(k)F(q)\q) L \k) R ® 6L fe 6+|0> 6 + F(k)F(q)\k) h \q) R ® bl q b{\0) b 



(Al) 



with ATlr, being a normalization constant, 
Nt R = (G^\V LR \GW) = (G 2 (fc) + G 2 (q)) (l + ±S\ + {\F(k)\ 2 + |F(g)| 2 ) 2 - 2G 2 (k)\F(k)f - 2G 2 (q)\F(q)f . (A2) 



We already introduced the sum S as being 

S = J2G 2 (k>). (A3) 

k' 

Appendix B: Continuum limit 

When taking the sum over all two-dimensional in-planc 
wave vectors k in Eqn. ( |44[ ), the appearing vectors de- 
pend on the boundary conditions. We choose peridodic 



boundary conditions, and hence 

k x = —n Xl n x = 0,±1, 



2n 



fly j fly 0, ibl , . . . , 



(Bl) 
(B2) 



where L x r y \ is the cavity length in a;(?/)-direction and 
n x ( y ) is an integer. Every discrete wave vector k has a 
volume A in fc-space: 



A = Ak x Ak y 



(27T) 2 = (2nl 

L X Ly 



A 



(B3) 
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and A.k x r y ) is the difference between two adjacent wave 
vectors in x(?/)-direction, A the sampie area. In the con- 
tinuum limit, the k vectors lie close in the reciprocal 
space and the sum can be replaced by an integral 

S = £ G\k) = I £ AG 2 (fc) -> I / d 2 fc G 2 (fc) 

7„ 7„ w 



2tt 



dkkG 2 (k) 



k=0 



27T V W 12, 



(B4) 



where we use polar coordinates to evaluate the integral 
and carried out the polar-angle integration. In the last 
step, we make a substitution and introduce the dimen- 



sionless variable k := — k to get the dimensionless inte- 

W12 ° 

gral 1, 



1= / dfcfcG 2 ( — k 

J \ W 12 
k=0 



(B5) 



One can show that the expansion coefficient G(k) de- 
creases like 1/fc 2 for large k. Consequently, the integrand 
has the asymptotics 



kG 2 {k) 



k— >oo 1 



(B6) 



and hence, the integral converges. We evaluate I nu- 
merically using the same parameters as above and find 
1= 1.9 -10 -4 . 
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